Inhomogeneity effects in oxygen doped HgBa2Cu04 
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We theoretically investigate inhomogeneity effects on the charges, electric field gradients and site- 
projected densities of states in HgBa2Cu04+,5. We find pronounced differences in the doping-induced 
number of holes at different atomic sites. The contributions of these sites to the density of states 
in the vicinity of the Fermi level are peaked at the same energy, but vary in magnitude by up to 
70 percent and have different energy dependence. Due to this energy dependence the role of the 
intrinsic inhomogeneities for superconductivity strongly depends on the energy and character of the 
quasiparticle mediating the Cooper pairing. Our results can explain the origin of doping-induced 
effects observed either by local or macroscopic experimental probes. 



PACS numbers: 74.81.-g, 74.25.Jb, 74.72.-h 

I. INTRODUCTION 

Doping is the most crucial parameter to influence the 
critical temperature in high T c cuprates. The metallic 
state is reached by doping through the replacement of 
ions or by introduction of excess oxygen. Thereby the 
dopant adds states at the Fermi level Ep, where these 
new charge carriers (holes) cause a downward shift of 
Ep and the carrier concentration in the Cu02 plane is 
changed. 

Most of the theoretical models make use of two simpli- 
fications to describe doping effects in the complex physics 
of these systems: First, they use the carrier concentration 
rather than the doping level as an input parameter, and 
second, they assume that the additional charge is uni- 
formly distributed within the different sites of the same 
atomic species in the Cu02 plane. Concerning the former 
point, it has been showni for HgBa2Cu04+<5 how doping 
influences the charge carrier distribution and what limits 
the amount of holes that can be created. For the lat- 
ter point, no first-principles calculations have been per- 
formed for any of the cuprates. 

At the same time there is strong experimental evidence 
for intrinsic inhomogeneous charge distribution in high 
T c compounds. For example, the existence of stripes 2 
in underdoped La2Cu04 and Bi based compounds is a 
clear manifestation of this. However, there are other 
sources of spatially non-uniform charge density. Scanning 
tunneling microscopy (STM) experiments performed on 
Bi2Sr2CaCu208+x exhibit an inhomogeneous surface car- 
rier distribution on a length scale of 14 A, in the normal 
as well as in the superconducting state, 3 It was claimed 
that this inhomogeneity has an unknown origin but is 
an intrinsic property of this doped material, i.e. not re- 
lated to impurities. More recently, it was found that the 
nanoscale electronic disorder in cuprates can be traced 
back to inhomogeneities on the atomic scale. 4 Local bulk 
probes like nuclear quadrupole resonance (NQR) and nu- 
clear magnetic resonance (NMR) experiments exhibit not 



only changes in the hole content with doping, but also the 
existence of different sites of the same species as it was 
seen in Tl based compounds^ and Sr doped I^CuO^ 
Since the spatial scale of inhomogeneities is determined 
by the amount of dopants, for doping levels in the regime 
of underdoped to optimally doped, investigations on a 
scale of a few lattice parameters are required. Wang 
et al. considered the role of the Coulomb impurities 
located close to the CuC>2 plane on the spatial depen- 
dence of the gap in a <i-wave superconductor within the 
t — J modelA In this work, we perform parameter-free 
calculations which allow a detailed analysis of the charge 
distribution. By this way the inhomogeneities can be 
quantitatively studied and linked to experimental find- 
ings. We investigate the doping-induced charge redistri- 
bution in HgBa2Cu04+5, which exhibits the highest T c 
among all single-layer cuprates. Focusing on the Cu02 
planes, we will show that the inhomogeneity is quite pro- 
nounced, it occurs on a scale of a few lattice constants, 
and is strongly doping dependent. We will also consider 
the consequences for the site-projected densities of states 
and their impact on superconductivity. To this extent, 
we have performed a series of supercell calculations cor- 
responding to oxygen concentrations of 5=1/8, 1/6, 2/9, 
and 1/4. 



II. SUPERCELL CALCULATIONS 

Unit cells of the 8-, 6-, 9-, and 4- fold size compared to 
the undoped case are presented in Fig. ^ The partic- 
ular environment with respect to the dopant makes the 
various copper and oxygen sites within the Cu02 planes 
inequivalent. The coordinates of these different types of 
Cu and O atoms are provided in Table [I] for the struc- 
tures shown in the two upper rows of the figure. The 
5 = 0.167 case gives rise to an orthorhombic cell with 
3x2 single unit cell volumes, the three other structures 
are tetragonal. 
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The charge of a given site calculated within the su- 
percell approach results from the influence of all dopant 
ions which may be generally long-ranged. However, due 
to the two-dimensional character of the carriers in the 
high-T c cuprates, the interaction of the dopant with its 
surrounding is screened at an in-plane length scale which 
is approximately the distance between the dopant and 
the Cu02 plane. 8 Hence this interaction is short-ranged. 
Therefore the sizes of the chosen supercells are big enough 
to study the effects of charge redistribution caused even 
by a single dopant ion. The fact that the doping-induced 
charges averaged over the respective unit cell are only 
very weakly sensitive to the choice of the supercell has 
already been pointed out in Ref. GJ In order to get more 
insight into the effects on the local environment, we con- 
sider an alternative supercell for 5 = 0.25 as given in the 
lowest panel of Fig. [I] This structure exhibits a less 
homogeneous distribution of dopants. Since it, however, 
is energetically unfavorable, as will be discussed below, 
the corresponding structural data and results are not dis- 
played in Table [H 

All calculations have been carried out within den- 
sity functional theory employing the full-potential lin- 
earized augmented plane- wave (LAPW) method utilizing 
the WIEN2k code. 9 Exchange and correlation effects are 
treated within the local density approximation (LDA). 
The computational details are given elsewhere. 1 



III. RESULTS 

A. Charge distribution 

The atomic-like basis functions used within the LAPW 
method inside the so-called atomic spheres 10 allow to an- 
alyze the charge density within each sphere by decom- 
posing it with respect to the contributions of different 
orbital characters. Although these numbers depend on 
the choice of the atomic-sphere radii they provide useful 
information about the charges around different atomic 
sites. 

The decrease in the partial charges in the atomic 
spheres averaged over the respective unit cells have been 
described in detail in Ref. 1. It was observed that a lin- 
ear increase in the total number of holes for 5 < 0.22, 
is followed by a plateau at higher 5. The reason for this 
saturation was found in the fact that at 5 ~ 0.22 the 
dopant reaches a closed shell O 2- configuration. The ra- 
tio of the additional average hole contents on Cu(d x 2_ y 2) 
and 0(p x ) orbitals hcu/ho is approximately 4.0. (Note 
that there are two oxygen atoms in the Cu02 plane per 
chemical formula unit, therefore the amount of oxygen 
holes is half that of copper.) 

Now we focus on the charge distribution on the atomic 
scale. Table U presents the number of holes h s created at 
the specific copper and oxygen sites in the Cu02 planes 
with the undoped material taken as the reference. As 
one can see, h s strongly depends on the atomic species, 
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FIG. 1: (Color online) Copper-oxygen planes of the 
HgBa2Cu04+<5 supercells representing the doping concentra- 
tions of 5=0.125 (upper left right panel), 0.167 (upper right 
panel), 0.22 (lower left panel), and 0.25 (lower right panel). 
The unit cells are indicated by the bold gray lines, the cop- 
per (oxygen) atoms by the red (blue) spheres. The projection 
of the dopant oxygen (located at the Hg plane) onto these 
planes is marked by the gray circles. In the lowest panel, an 
alternative structure for 5=0.25 is displayed. 



the specific site, and the doping content. Concerning the 
doping dependence, for all given 6 values lower than 0.25, 
the maximal difference in the number of doping-induced 
holes among non-equivalent sites of the same species (de- 
noted as A/^ ax ) is approximately 0.01 e, although the 
distances from the ions to the next dopant may be quite 
different at smaller doping concentrations. For a fixed 
doping content, we define a measure for the inhomogene- 
ity, A/^ ax divided by the average doping-induced charge 
of this species, i.e. hc u or ho. This value decreases for 
copper from 21% for 5=0.125 to 13% for 5=0.22 (and zero 
for 5=0.25). For oxygen it is much more pronounced with 
values of 114%, 95%, 51%, and 14% for 5=0.125, 0.167, 
0.22 and 0.25, respectively. Summarizing these findings 
we conclude that the relative inhomogeneity with respect 
to the total amount of doping-induced holes (i) is more 
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TABLE I: Copper and oxygen types, their multiplicities, their 
(X,Y) coordinates given in lattice constants A and B of the re- 
spective supercells, and their in-plane distances to the nearest 
dopant oxygen do, where a is the single cell lattice parameter. 
h s denotes the number of holes created by doping at the given 
site, V 11 are the diagonal components of the electric field gra- 
dient tensor in 10 21 V/m 2 units, and r\ is the corresponding 
asymmetry parameter. The principal axes of the EFG are 
either identical or close to the crystal axes. 



pronounced for oxygen than for copper, since doping in- 
troduces less holes on the oxygen than on the copper 
sites, and (ii) is decreasing with increasing S. 

The isotropy of the charge distribution for S = 0.25 
given in the table may be considered as an artifact due to 
the choice of the supercell. For this reason, we have per- 
formed calculations for an alternative structure, which is 
displayed in the lowest panel of Fig. [I] For this case, the 
inhomogeneity is more pronounced, with Aa™ ax = 0.013 
for copper, still being in the range of a hundredth of an 
electron as found for the other doping concentrations. 
But more important, the total energy of this structure is 
more than 5 mRy per formula unit higher, which corre- 
sponds to a temperature of about 900K. This large energy 
compared to the other cell choice can be explained by the 
smaller distance between the nearest neighbor dopants, 



which, being highly negatively charged, repel each other. 
For the other doping concentrations, we have chosen the 
supercells such that the separation of dopants is maxi- 
mal. 

The non-uniform charge distribution can also be un- 
derstood in terms of the distances to the next dopant 
ion. Focusing on the projection onto the Cu02 plane the 
spacings to the copper positions are ay/2/2, avTO/2, etc. 
Although for all 5 values the maximum amount of holes 
at Cu is created in the first coordination shell with the 
distance being ay/2/ 2^ its value strongly depends on S. 
While for S = 0.125 and 0.167 there are only two non- 
equivalent copper positions differing in the distance to 
the dopant, there are four non-equivalent Cu atoms for 
S = 0.22. In this case, the biggest amount of holes is 
created at Cul, which is influenced by two dopant sites 
with the same distance; Cu2 exhibits the same spacing 
to the nearest dopant atom, but the next nearest dopant 
oxygen is farther away. A similar analysis can be done 
for the oxygen sites, where the in-plane distances to the 
dopant are a/2, ay/5/2, 3a/2, etc. For example, at S = 
0.125, 03 (at a distance of a/2) is most affected by the 
excess oxygen, followed by 02, which interacts with a 
dopant at a distance of ay/5/2, and Ol where the dopant 
is 3a/2 away. Also for the oxygen sites the amount of 
created holes very strongly depends on the oxygen con- 
centration, when positions with the same spacing to the 
dopant are compared. The values for the nearest dopant 
range from 0.0139 for S = 0.125 via 0.0195 for S = 0.167 
to 0.0241 for S = 0.22. 



B. Electric field gradients 

The distribution of doping-induced charges can be 
probed by nuclear quadrupole and magnetic resonance. 
The electronic and nuclear charge densities form a non- 
uniform electric field within the unit cell. A nuclear 
quadrupole moment is coupled to the gradients of the 
electric field 11 given by V a p = dE a /dxp, where E a is 
the component of the local electric field, and xp is the 
Cartesian coordinate. The traceless tensor V a p can be di- 
agonalized and written in the principal axes as V 11 . The 
isotopes 17 O and 63 Cu are usually used as the probes of 
electric field gradients in high-T c cuprates. As a result 
of the interaction of their nuclear quadrupole momenta 
Q with the electric field gradients at the lattice sites, a 
resonance line at 

^=^ z (l + r ? 2 /3) 1/2 (3.1) 

appears in the radio frequency wave absorption spectrum 
of the system. Ordering the tensor elements by their ab- 
solute value, such that V xx is the smallest and V zz is 
the maximal component, the asymmetry parameter 77 is 
defined as (V xx -V yy )/V zz , where V zz is called the elec- 
tric field gradient EFG. The corresponding tensor com- 
ponents and asymmetry parameters are given in Table 
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[IJ Some of the site symmetries give rise to non-diagonal 
components, which are, however, very small. Therefore 
the principle tensor axes nearly coincide with the crystal 
axes and thus V 11 are given as the corresponding diag- 
onal elements. Recent experimental data 6 on EFGs in 
La2-^Sr :c Cu04 clearly show the following observations: 
There is (i) a doping-induced upward shift of the mean 
copper EFG in the order of 10%, (ii) a splitting of the res- 
onance line according to different inequivalent sites, and 
(iii) a decrease in the splitting of the resonant frequencies 
of two different Cu sites with increasing doping. As one 
can conclude from Table [IJ all these findings are consis- 
tent with our results™ NMR measurements on Tl-based 
compounds 5 exhibit changes in the Cu and O EFGs of 
the order of 10% with doping, again concomitant with 
our calculations. At the same time, NQR data taken 
on Hgl202 also revealed an upward shift with doping, 13 
which is, however, too large in order to be understood in 
terms of oxygen-doping induced holes only. We suppose 
additional defect s 14 ! 15 to be responsible for such a strong 
increase of the resonant frequencies. Our results are in a 
good qualitative agreement with the experimental obser- 
vations of Ref. [l6l performed on under- and overdoped 
HgBa2Cu04 +( 5 crystals. In these measurements, a shift, 
a splitting, and a change of the in-plane 17 O NMR line 
width were revealed with the change of the doping con- 
centration. 



C. Density of states 

Having discussed the site-dependent properties, we 
turn our attention to the density of states (DOS), which 
is an integral characteristic of the system and should be 
crucial for the transition temperature. In Fig. the site- 
projected densities of states are presented for £=0.167 
and 0.22, highlighting the contributions of selected cop- 
per and oxygen spheres. Like the charge carriers, the 
DOS's also exhibit pronounced inhomogeneities which 
are of similar order of magnitude. The relative differ- 
ences in the DOS at the Fermi level from different oxygen 
sites are up to roughly 70 % in both cases, while these 
are much smaller for the copper spheres. For £=0.167 
the two values are nearly the same, while the differences 
are somewhat bigger for £=0.22, where the smallest and 
largest contribution are 0.468 (not shown in the figure) 
and 0.500, respectively. This feature is much more pro- 
nounced when moving away from the Fermi level in the 
range of 0.05 eV which is a typical energy scale of quasi- 
particles mediating the pairing. At (5=0.167 all contri- 
butions show broad peaks at the same energy which is 
slightly below the Fermi level. For (5=0.22 the peaks have 
become much sharper and have moved somewhat above 
Ef. Again, all copper and oxygen contributions exhibit 
their maxima at the same position. 

Although charges and densities of states demonstrate 
similarities, which are inhomogeneities of the same order 
of magnitude, there is a clear difference in their behavior 
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FIG. 2: (color online) Selected Cu(d x 2_ y 2) and 0(p x ) con- 
tributions to the density of states in states per eV stemming 
from the respective atomic spheres for two different doping 
concentrations. The vertical lines indicate the Fermi level. 



when we focus on the DOS in the vicinity of the Fermi 
level. This comes from the fact that the charges are ob- 
tained by integrals of the DOS over the whole valence 
energy range. Yet, only a small energy range around Ep 
is relevant for superconductivity. A noticeable fact in 
this context is that upon doping the pda band becomes 
less populated, which is important for the charges in the 
Cu02 plane. Regarding the DOS at Ef-, also the posi- 
tion of the saddle point (often referred to as van Hove 
singularity) with respect to the Fermi level plays a role. 



IV. CONCLUSIONS 

To summarize, we have performed ab-initio calcula- 
tions for the oxygen doped HgBa2Cu04 to consider in- 
homogeneity effects on the charges and densities of states. 
We find pronounced differences in the doping-induced 
number of holes at given copper and oxygen sites, where 
the degree of the inhomogeneity strongly depends on the 
doping level and the atomic species. Our results are con- 
sistent with the findings of Pan et al& giving a spatial 
scale of a few lattice parameters. Furthermore, they ex- 
plain the behavior of NMR and NQR lines, i.e. a shift 
and a splitting with doping^Sii^ 

Besides the atomic-scale experiments discussed above, 
there could be further consequences of the non-uniform 
charge distribution on larger spatial scales. This distri- 
bution gives rise to a random potential for the carriers 
near the Fermi level leading to their momentum relax- 
ation. On the macroscopic scale, these effects have been 
seen in the resistivity measurements^ and electronic Ra- 
man spectroscopy™ In a d x i_ y i superconductor, such 
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a potential is pair-breaking and leads to a reduction of 
T c . This decrease, however, turned out to be smaller 
than expected from the momentum relaxation rate, as 
it was shown in Ref. On the nanoscale, this ran- 
dom potential leads to the electron inhomogeneity for 
the low-energy states probed by the scanning tunneling 
spectroscopy in the superconducting state A2fii2ii2i 

Concerning the density of states, which is directly con- 
nected to superconductivity, we observe that the contri- 
butions of all plane copper and oxygen species are dif- 
ferent in magnitude but peak at the same energy. Going 
below and above Ep we find more pronounced differences 



which makes us conclude that the importance of the in- 
homogeneities for superconductivity strongly depends on 
the energy and character of the quasiparticle mediating 
the Cooper pairing. 
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